The Effect of Large Amplitude Fluctuations in the Ginzburg-Landau Phase 

Transition 
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The lattice Ginzburg-Landau model in d = 3 and d = 2 
is simulated , for different values of the coherence length £ in 
units of the lattice spacing a, using a Monte Carlo method. 
The energy, specific heat, vortex density v, helicity modulus 
T M and mean square amplitude are measured to map the 
phase diagram on the plane T — £. When amplitude fluctua- 
tions, controlled by the parameter £, become large (£ ~ 1) a 
proliferation of vortex excitations occurs changing the phase 
transition from continuous to first order. 

PACS numbers: 74.25.Dw, 64.60.-i, 05.70.Fh 

The Ginzburg-Landau (G-L) model, involving a 
complex field ip — \ip\ exp[i6], is a kind of multi purpose' 
or metamodel' which captures the essence of several in- 
teresting phenomena of condensed matter: superfluid- 
ity, superconductivity and the melting transition. It 
turns out that in all these phase transitions vortex-like 
defects (topological singular phase configurations) play 
a central role condensing above the critical temperature 
T c and breaking the ordered state existent below T c JjJ] . 
With the discovery of high-temperature superconduc- 
tors an entire new area in condensed matter devoted 
to vortex physics has been opened Q , [gj . A very rich 
phase diagram for the vortex matter emerged both from 
the experimental and theoretical side. 

The nature of the G-L phase transition in d=3 
and d=2 systems still remains controversial theoreti- 
cally. For d=3, in the early 70's, Wilson and Fisher 
using renormalization group techniques, concluded that 
the G-L model belongs to the same universality class of 
the XY model ("phase only" approximation), which is 
known to exhibit a second-order phase transition. On 
the other hand, at T — T c amplitude fluctuations, con- 
trolled by the coherence length in units of the lattice 
spacing £, are in principle not negligible and they might 
affect the critical behavior Q. In that sense, a recent 
variational approach |(| showed evidences that for d =3 
the G-L transition can turn first order due to the inter- 
play between phase and amplitude fluctuations when 
the latter are large enough. For d — 2 doubt remained 
if the first order found was not an artifact of the used 
approximation. However, for d = 2, we found evidences 
of a first order transition when amplitude fluctuations 
are allowed to be large by a suitable choice of the G-L 
hamiltonian parameters [J7j. 

In this letter, using a Monte Carlo method we show 



that: 1) the G-L phase transition becomes first order 
when £ is greater than a value £2^1 close to 1, allowing 
non-negligible amplitude fluctuations and 2) that this 
change of order is connected with a sudden prolifera- 
tion of vortices. We worked on square lattices of size L 
and spacing a, we denote the lattice sites by x and the 
lattice links by (x, fx) with fx = 1, .., d. We used the pa- 
rameterization of the Boltzmann exponent j3H of ref. J5j 
i.e. in terms of two parameters, the coherence length £ 
and a lattice temperature Tl, and a dimensionless order 
parameter -0. This parameterization is connected with 
the ordinary Wilson parametrizarion H\y /T H by: 
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where 



is the constant value to what 



ijj approaches infinitely deep in the interior of the super- 
conductor. In what follows we will omit the subscript 
L of T. In the limit of £ = (or u = 00) the radial 
degree of freedom is frozen and the model reduces to 
the XY model. The limit we are interested in is just 
the opposite: £ ~ 1, where large amplitude fluctuations 
occur. 

The calculations were performed using periodic 
boundary conditions (PBC). As it is common practice, 
we have discretized the 0(2) global symmetry group to 
a Z(N) to increase the speed of the simulation. For 
Z(60), we found no appreciable differences when com- 
paring results with previous runs carried out with the 
full 0(2) group in relatively small lattices. We used a 
sequence of lattice sizes from L=10 to L=64 for d = 2 
and from L—6 to L—16 for d — 3. For the case up 
to N = L d = 508 sites, we thermalized with usually 
20.000-40.000 sweeps and averaged over another 50.000- 
100.000 sweeps. For greater N larger runs were per- 
formed, typically 50.000 sweeps were discarded for equi- 
libration and averaged over 150.000-250.000 sweeps. We 
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wish to emphasize that for the case £ ~ 1, we are most 
interested in, no appreciable differences take place by 
taking N > 508 (for instance 8 3 and 16 3 give very sim- 
ilar results). The errors for the measured observables 
are estimated in a standard way by dividing measures in 
bins large enough to regard them as uncorrelated sam- 
ples. The following observables were measured to map 
the phase diagram and analyze the phase transition: 

i) The energy density e =< H > /N and the spe- 
cific heat cy. cy was computed both simply as the en- 
ergy variance per site i.e. cy = (< H 2 > — < H > 2 )/N 
and as the temperature derivative of e. 

ii) The vortex density v (density of loop vortices 
in d=3 and point vortices in <i=2), which serves as a 
measure of the phase disorder, defined by: 

v = j2^2\ m *p\' ( 3 ) 

where *p denotes the lattice cell dual to a given lattice 
plaquette p (a site in d = 2 and a link in d = 3) and the 
quantity m* p is the vortex "charge" assigned to *p and 
measured over the plaquette p as: 

m* P = - B 2 ] 2lT + [0 2 - 8 3 } 27r + 

+ [03-04] 27r + [04-0l] 27r ), (4) 

where [a]27r stands for a modulo 2n: [ch^tt — ol + 2im, 
with n an integer such that a + 2im G (— 7T,7r], hence 
™* P = n>i2 + n23 + n.34 + 7141 can take three values: 0, 
±1 (the value ±2 has a negligible probability). 

iii) The helicity modulus T M [|| T M measures the 
phase-stiffness along the direction \i. For a spin system 
with PBC the helicity modulus measures the cost in free 
energy of imposing a "twist" equal to L5 in the phase 
between two opposite boundaries of the system. T M is 
computed using the expression Jl0| : 

X 

~Y K E I ^ I' ^*+" I sin (<W - e *)f >}■ ( 5 ) 

x 

iv) The mean square amplitude p =< > which 
takes its minimum value at the phase transition. 

First, in Fig. 1 we report the phase structure of 
the model in the (T, £) plane for d=3 (o) and d—2 (•). 
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FIG. 1. Phase diagram in the (T, £) plane for d—2 (•, 
L = 40) and d=3 (o, L = 12). 



For £ larger than £2^1 the order of the phase tran- 
sition changes from continuous (dashed line) to first- 
order (filled line). Specifically, £2^1 — 0.8 {d=2) and 
£2^1 — 0.7 (d=3). For £ > £2—1 £ exhibits a large 
hysteresis phenomenon and an histogram with a dou- 
ble peak structure, corresponding to the two coexisting 
phases, both features characteristic of a first-order tran- 
sition. In Fig. 2 we present energy histograms for d=2 
for 2 values of £: £=0.85 > £2^1 and £=0.75 < £ 2 ->i, 
showing the different behavior as L increases: both 
peaks remain fixed as L increases in the first case while 
they approach each other as L increases in the second 
case (besides the peaks are lower and wider). Further- 
more, for £ = 0.85 the width of each of the peaks clearly 
scales as ^/[ttt) = \ ■> due to ordinary non-critical fluc- 
tuations, while for £ = 0.75 the peaks do not scale this 
way. (Similar histograms are found for d—3 for £ > £2^1 
and £ < £2^1). 
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FIG. 2. Histograms of e in d=2 for L=10 (below) ,L=20 
(middle) and L=40 (above), (a) £=0.85: the two peaks be- 
come sharper and their position remain fixed as L increases, 
(b) Zoom of the right peak showing the scale of its width as 
1/L. (c) £=0.75: the width of the two peaks do not scales 
as 1/L and they approach each other as L increases. 



Fig. 3 shows the hysteresis phenomenon found when 
considering heating and cooling runs for d=3 (a com- 
pletely similar hysteresis diagram is found for d=2 for 
£>^i). 
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FIG. 3. Hysteresis in d=3, L=14 for s, £ = 1 (error bars 
are smaller than the symbol size) . 



In order to illustrate the central role played by vor- 
tex excitations in triggering and determining the nature 
of the phase transition, Fig. 4-(a) and 4-(b) show a plot 
of f vs. T for different values of £ respectively for d=2 
and d=3. For coherence lengths larger than £2-»i we ob- 
serve a sharp jump in the vortex density v evidencing 
a sudden proliferation of vortices which coincides with 
the discontinuity in s. As long as we decrease £ the 
jump becomes more smooth and moves to higher val- 
ues of T c . In particular, for d=2, for £ ~ 0.1 something 
very close to the Kosterlitz-Thouless (K-T) behavior is 
reached. The increase in the density of vortices when 
amplitude fluctuations are large is due basically to the 
fact that amplitude fluctuations decrease the energy of 
vortices enhancing vortex production. The same hap- 
pens for the XY model with modified interaction O]; 
in fact, the shape modification of the interaction can be 
connected with a core energy variation. 
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FIG. 4. Plot of v vs. T (error bars are smaller than the 
symbol size. Left: d=2, L=40; £=1 (A), £=0.8 (V), £ = 0.5 
(x), £ = 0.1 (+) and XY (o). Right: d=3, L=12; £=1 (A), 
£=0.8 (V), £=0.7 (o) and £=0.5 (x) 
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FIG. 5. Plot of < |^| 2 > vs. T. (a) d = 2: £=2 ([ 
£=1 (A), £=0.8 (V), £=0.7 (*) and £=0.5 (x). (b) d = 
£=2 (□), £=1 (A), £=0.8 (V), £=0.7 (*) and £=0.5 (x). 
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Figures 5 shows a plot of p =< \ip\ 2 
different values of £. Note that for £ > £2 
a sharp drop which coincides with the jump in v and e 
(and the drop in not depicted). 

Therefore, the nature of the phase transition of the 
G-L model depends dramatically on the value of the 
coherence length £. For £ > £2^1 the transition to a 
disordered state implies latent heat, a a discontinuous 
jump of the vortex density and the subsequent abrupt 
drops in and p i.e. all the features of a first order 
transition. On the other hand, for ( « 1 the G-L re- 
duces to the XY model (with the more subtle K-T phase 
transition for d=2). For d=3, the value of £2^1 — 0.7 is 
in agreement with ref. |J in which the authors found the 
ff-^i > 1/4.5 For d=2, the value of £2—1 — 0.8 belongs 
to the region (£ ^1) found im ref. |J such that the RG 
trajectories cross the first order line of Minnhagen's |0] 
generic phase diagram for the two-dimensional Coulomb 
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gas. 

To conclude, both for d=2 and d=3, we offer clear 
evidences that the order of the phase transition, trig- 
gered by vortices, in the G-L model depends on the 
value of £ which controlls the size of amplitude fluctua- 
tions. Whether or not such mechanism based on the in- 
terplay between amplitude and phase fluctuations takes 
place in real systems like high T c superconductors or in 
explaining the change in order of the melting transition 
is something which deserves a more careful analysis. 
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